Influence of light nuclei on neutrino-driven supernova outflows 
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We study the composition of the outer layers of a protoneutron star and show that light nuclei are 
present in substantial amounts. The composition is dominated by nucleons, deuterons, tritons and 
alpha particles; ^He is present in smaller amounts. This composition can be studied in laboratory 
experiments with new neutron-rich radioactive beams that can reproduce similar densities and 
temperatures. After including the corresponding neutrino interactions, we demonstrate that light 
nuclei have a small impact on the average energy of the emitted electron neutrinos, but are significant 
for the average energy of antineutrinos. During the early post-explosion phase, the average energy 
of electron antineutrinos is slightly increased, while at later times during the protoneutron star 
cooling it is reduced by about 1 MeV. The consequences of these changes for nucleosynthesis in 
neutrino-driven supernova outflows are discussed. 
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I. INTRODUCTION 

Core-collapse supernovae are the birth places of neu- 
tron stars. During the explosion, around 10^'^ ergs, cor- 
responding to the binding energy of the neutron star, are 
emitted in all neutrino species from a thermal surface de- 
noted as neutrinosphere. The spectrum and luminosities 
of the neutrinos radiated from the newly formed neu- 
tron star, the so-called protoneutron star, are of great 
importance for the interpretation of supernova neutrino 
detections, for studies of flavor conversion in the star 
mantle and for nucleosynthesis occurring in the 

so-called neutrino-driven wind [3, Q ■ 

A crucial input for protoneutron star evolution is the 
equation of state (EOS). Currently available simulations 
use either the Lattimer and Swesty Q or Shen et al. Q 
EOS, or an EOS based on nuclear statistical equilibrium 
(NSE) f?\. All these EOS describe the nucleonic compo- 
sition by a mixture of neutrons, protons, alpha particles 
and a representative heavy nucleus. While these EOS 
include the essential composition for the relatively low 
densities and high temperatures present in the protoneu- 
tron star atmosphere, where matter is fully dissociated, 
they cannot account for the composition in the region 
of the crust where the larger densities [p ~ lO^^g/cm^) 
allow for the formation of light nuclei with A = 2 and 3 
in addition to alpha particles [1, Q . In contrast to the 
conditions in cold neutron stars, the high temperatures 
(T > 4 MeV) of the protoneutron star crust suppress the 
formation of heavier nuclei. Deeper in the protoneutron 
star the densities become so large that light nuclei melt 
and a transition from inhomogeneous phases to homo- 
geneous nuclear matter takes place [HI • As we will 
discuss below, the change in composition caused by the 
presence of light nuclei in the outer crust affects the neu- 



trino opacities and consequently changes their spectra 
and luminosities. In addition, light nuclei are present in 
the region behind the shock, where the emitted neutrinos 
are expected to deposit their eiiergy in the delayed super- 
nova explosion mechanism [T^ . Il3j | and therefore interac- 
tions with neutrinos have to be included as well [1, [13] . 

The neutrinos emitted by the young protoneutron star 
produce an outflow of baryonic matter known as the 
neutrino-driven wind that has been the subject of many 
studies including full hydrodynamical simulations [7| , an- 
alytical [1| and steady-state approaches (see Ref. and 
references therein). This outflow is initially very hot and 
essentially consists of free neutrons and protons in a ra- 
tio that is determined by the competition of neutrino and 
antineutrino absorptions on nucleons and their inverse re- 
actions. But as the matter expands and cools, nucleons 
can be assembled into nuclei, and elements even heavier 
than iron can be formed. If this occurs with a large abun- 
dance of free neutrons present, these can be captured on 
heavy nuclei (the "seed" ) producing an r-process [l^, [l3] ■ 
For a successful r-process a large neutron-to-seed ratio is 
necessary, requiring outflows with short dynamical time 
scales (a few milliseconds), high entropies (above 150 fee) 
and low electron fractions (Ye < 0.5) 18]. 

In the present paper, we explore the influence of light 
nuclei on the spectra and luminosities of electron neutri- 
nos and antineutrinos emitted during the cooling phase of 
the protoneutron star and consequently on the electron- 
to-baryon ratio of the ejected matter. The existence of 
light nuclei potentially also affects ji and r neutrinos. We 
plan to investigate this in future work. The early post- 
bounce evolution and pre-explosion phase might also be 
affected by the presence of light nuclei in the matter com- 
position. However, this issue cannot be explored here 
in detail, because the neutrino transport conditions in 
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those phases and the development of hydrodynamical in- 
stabiHties in the forming neutron star and in the neutrino 
heating region (for reviews, see Refs. (H,!!^) require full 
radiation hydrodynamics simulations, which are beyond 
our study here. 

If light nuclei are present in the region near the ncutri- 
nospheres they will influence neutrinos and antincutrinos 
differently. Deep in the interior of the protoneutron star 
neutrinos, are in chemical equilibrium with matter. At 
larger radii the decrease of temperature and density al- 
lows for neutrinos to decouple from matter near the pro- 
toneutron star surface and near the region where nuclei 
form. For the very neutron-rich conditions present in this 
environment, the formation of light nuclei occurs mostly 
at the expense of free protons. As protons represent the 
major source of opacity for electron antincutrinos, and 
because the antincutrino cross sections on protons and 
light nuclei are different for the relevant energies, light 
nuclei can potentially affect the spectra and luminosities 
of the emitted antincutrinos. The situation is not the 
same for electron neutrinos. The abundance of neutrons 
is so large that it is insignificantly changed by the ap- 
pearance of light nuclei. Consequently, the electron neu- 
trino opacity, which is dominated by interactions with 
neutrons, remains practically unchanged. This asymmet- 
rical effect on the radiated neutrinos and antincutrinos 
can potentially change the proton-to-neutron ratio of the 
ejecta and consequently the nucleosynthesis. A study of 
these effects requires an EOS that includes light nuclei in 
the composition and their corresponding neutrino cross 
sections. 

This warm nuclear system near the neutrinosphere, 
including the abundances of light nuclei, can be stud- 
ied with heavy- ion collisions. Fragments emitted from a 
system at intermediate velocities may come from a low- 
density region between the colliding nuclei. The density 
and temperature of this region can be similar to condi- 
tions near the neutrinosphere. Kowalski et al. have 
measured the abundances of light nuclei (deuterons, tri- 
tons, '^He, and ^He) in near Fermi energy heavy-ion col- 
lisions of ^^Zn on ^^Mo and on -^^''Au. They found ^He 
abundances and symmetry energies similar to those pre- 
dicted by the virial EOS [2l|. In the future, experiments 
with radioactive beams will enable studies of the more 
neutron-rich neutrinosphere conditions. 

This paper is organized as follows. Section |TT] 
presents the EOS (which we use for the determination of 
light-element abundances), the neutron star atmosphere 
model, and the neutrino cross sections with light nuclei. 
Section IIIII discusses the impact of light nuclei on the 
spectra and luminosities of electron neutrinos and an- 
tincutrinos and their influence on the electron fraction of 
the ejecta. Finally, we conclude in Sect. II VI 
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FIG. 1: Neutron star atmosphere profiles of density and tem- 
perature corresponding to model M15-ll-rl of Ref. 01 for 
times t = 2,5,7 and 10 s post bounce. 



II. EQUATION OF STATE, OPACITIES, AND 
NEUTRINOSPHERE DETERMINATION 

A. Neutron star model and equation of state 

Our work employs the protoneutron star model M15- 
11-rl of Ref. 0. The model describes the structure of 
the surface layers and of the neutrino-driven wind of 
a protoneutron star with a baryonic mass of 1.4 Mq, 
obtained in a spherically symmetric simulation of the 
(parametrized) neutrino-driven explosion of a 15 Mq 
star. The thermodynamical state of hot, dense mat- 
ter including its baryonic composition is fully character- 
ized by three independent variables, for example, density, 
temperature, and electron chemical potential (or baryon 
density, energy density, and electron fraction). In the su- 
pernova context, these have to be determined by solving 
the equations of hydrodynamics for the stellar plasma 
and the transport equations for neutrinos and antineu- 
trinos of all flavors, making use of an EOS that relates 
the pressure and all other thermodynamic variables to 
the three basic ones. Figure [T] shows the time evolution 
of the density and temperature profiles after core bounce 
in the region around the neutrinospheres ("neutron star 
atmosphere") as predicted by model M15-ll-rl. We point 
out that this model did not include light elements, there- 
fore the temperature and density are not exactly those 
one would obtain if light elements were included. 

Using these temperature and density profiles, we can 
determine the electron chemical potential (or, equiva- 
lently, the electron fraction Yg) by making the assump- 
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tion that the matter is in (neutrinoless) beta equihbrium. 
This is fairly well fulfilled in the neutron star atmosphere 
at sufhciently late post-bounce times when the delep- 
tonization of the protoneutron star interior has slowed 
down and the nascent neutron star radiates electron neu- 
trinos and antineutrinos with very similar number lumi- 
nosities. In this case the chemical potentials of neutrons, 
protons, and electrons fulfill the equality /i„ = fj,p + fj,e, 
which allows one to compute by invoking also charge 
neutrality, — Yp, and the relation for the total baryon 
number, Yp — 1~ Y^. Here Yp and Y^ denote the number 
fractions of free plus bound protons and neutrons, respec- 
tively. It is clear that the Ye-profile one thus obtains will 
depend on the considered EOS. 

We will compare results obtained with three differ- 
ent EOS. The first is the one used in Ref. 0], which 
essentially describes the baryonic composition as a non- 
interacting Boltzmann gas of neutrons, protons, and al- 
pha particles in NSE (a representative heavy nucleus is 
also included but of no relevance here). Calculations per- 
formed with this EOS serve as the reference to which we 
compare our results. The second EOS, denoted by NSE 
EOS, also assumes matter in NSE but consists of sev- 
eral thousand nuclei, for which partition functions have 
been computed in Ref. 241 , and includes Coulomb correc- 
tions [l^H^ and Fcrmi-Dirac statistics for neutrons and 
protons. Our third EOS is the virial EOS [1, [H, [H, [H , 
which is based on nuclei with ^ 4 and their inter- 
actions through second virial coefficients derived directly 
from scattering phase shifts. Figure [2] shows the composi- 
tion obtained from these different EOS for the thermody- 
namical conditions of Fig. [Tl assuming beta equilibrium 
for each of the displayed cases. In agreement with our 
arguments, we observe that in the vicinity of the neutri- 
nospheres the virial and NSE EOS predict deuteron and 
triton mass fractions significantly larger than the one of 
free protons. 

The virial and NSE EOS lead to very similar compo- 
sitions up to densities p ^ 10^'^g/cm^. At these densi- 
ties, the treatment of nuclear interactions in the virial 
EOS becomes unreliable, and in the NSE EOS interac- 
tions are neglected. This is signaled by a sudden increase 
in the abundance of heavy nuclei, when using the NSE 
EOS, or by a sudden drop of the proton mass fraction 
in the virial EOS, due to the breakdown of the virial ex- 
pansion with only second virial coefficients. For densities 
lower than p ^ 10^'^g/cm'^, we observe that the main 
differences of the virial and NSE EOS are in the alpha 
particle mass fractions due to attractive nucleon-alpha 
interactions |,21|] . 



B. Neutrino opacities and neutrinospheres 

To determine the neutrinosphere radius we follow 
Refs. [13, and define the effective neutrino opacity 
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FIG. 2: (Color online) Mass fractions (defined as the mass 
density of species i divided by the total mass density) for 
nuclei present around the surface of the protoneutron star 
using three different EOS: the one used in Ref. 7] (dashed 
lines^the NSE EOS (dotted lines), and the virial EOS 8, 2% 
[13, [13| (solid lines) . Beta equilibrium was assumed in each 
case. The lines labeled "heavy" represent the mass fraction 
of nuclei with ^4 > 4, which are included in the NSE EOS but 
not in the others. From top to bottom the profiles correspond 
to t = 2, 5, 7 and 10 s post bounce. The vertical lines mark the 
positions of the neutrinospheres of electron neutrinos (right 
line) and electron antineutrinos (left line). 
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for energy exchange or thermalization by 

KcS = V '«abs('«abs + ^scatt) • (1) 

The absorption opacity Kabs is considered to include all 
processes in which neutrinos exchange energy with the 
stellar medium, while the scattering opacity Kscatt con- 
tains those processes where mostly the momentum of 
the neutrinos is changed but essentially not their energy. 
These opacities, k = Y1 ntai, are obtained from the neu- 
trino cross sections ai and number densities n.i of the 
target particles in the stellar plasma. 

In order to determine an average neutrinosphere ra- 
dius, rather than an energy-dependent one, we assume 
neutrinos to be in thermal equilibrium with matter up to 
their so-called average energy sphere, where the bulk of 
the neutrino spectrum begins to decouple thermally from 
the background medium of the star. We will consider this 
energy sphere as the appropriate neutrinosphere in the 
context of the work presented in this paper. Until this 
location the neutrino phase-space distribution will be as- 
sumed to be a Fermi-Dirac distribution function for the 
local gas temperature. Outside of the neutrinosphere the 
spectral temperature is taken to be fixed to its value at 
the neutrinosphere. Thus, we average all opacities over 
the relevant neutrino spectral distributions and define: 

(Kcs) = (^abs) (('«abs} + (Kgcatt)) ■ (2) 

In the calculation of Kscatt we include elastic scattering 
off nucleons [2§] and nuclei [13] for both electron neutri- 
nos and antineutrinos. For determining Kabs for electron 
neutrinos, it is sufficient to consider their absorption on 
neutrons, because neutrons dominate the composition by 
far. In the case of electron antineutrinos, besides weak 
processes [2^ with the rare protons we have to include 
also inverse nucleon-nucleon bremsstrahlung ^sT] and in 
particular charged-current and neutral-current interac- 
tions with deuterons [s^l and tritons (as discussed be- 
low). The situation is different for '^He as target particle. 
Its abundance is always very low so that this contribu- 
tion to the antincutrino opacity can be neglected. More- 
over, inelastic scattering on electrons and positrons and 
neutrino-antineutrino annihilation can be neglected, be- 
cause either the cross sections or target densities of these 
reactions are typically smaller than those of the neutrino 
interactions with baryonic targets. Ignoring these pro- 
cesses has hardly any influence on the relative changes of 
the neutrinospheric positions that we intend to discuss 
in this paper. Our approach for estimating the influ- 
ence of composition effects in the surface layers of nascent 
neutron stars on the position of neutrino-matter decou- 
pling is rather qualitative and approximative anyway. 
It is certainly not suitable for making exact quantita- 
tive predictions, a goal that definitely requires radiation- 
hydrodynamics simulations with energy-dependent neu- 
trino transport. 

The first ab-initio calculations for neutral-current in- 
elastic cross section off '^H and ^He nuclei were pre- 
sented in Ref . Q , where the neutrino energy was averaged 



TABLE I: Cross sections for charged-current and neutral- 
current antineutrino scattering off "^H as a function of incident 
antineutrino energy E. 
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6.13 X 10"^ 
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9.32 X 10"^ 


1.17 X 10"^ 
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1.35 X 10"^ 
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1.90 X 10"^ 


2.05 X 10"^ 
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2.58 X 10"^ 


2.64 X 10"^ 
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3.43 X 10"^ 


3.34 X 10"^ 
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4.46 X 10"^ 


4.17 X 10"^ 
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5.70 X 10"^ 


5.14 X 10"^ 
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7.18 X 10"^ 


6.26 X 10"^ 
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8.93 X 10"^ 


7.54 X 10"^ 
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1.10 X 10° 


9.01 X 10"^ 
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1.61 X 10° 


1.25 X 10° 
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2.28 X 10° 


1.70 X 10° 
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3.13 X 10° 


2.24 X 10° 
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4.20 X 10° 


2.89 X 10° 
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5.50 X 10° 


3.66 X 10° 
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7.08 X 10° 


4.56 X 10° 
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8.97 X 10° 


5.60 X 10° 
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1.12 X 10^ 


6.78 X 10° 
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1.38 X 10^ 


8.13 X 10° 
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1.69 X 10^ 


9.63 X 10° 
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2.43 X 10^ 


1.31 X 10^ 
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3.38 X 10^ 


1.74 X 10^ 
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4.56 X 10^ 


2.23 X 10^ 
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5.99 X 10^ 


2.81 X 10^ 
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7.69 X 10^ 


3.45 X 10^ 
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9.65 X 10^ 


4.17 X 10^ 


78 


1.18 X 10^ 


4.96 X 10^ 
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1.42 X 10^ 


5.82 X 10^ 
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1.69 X 10^ 


6.75 X 10^ 


90 


1.98 X 10^ 


7.74 X 10^ 


95 


2.37 X 10^ 


9.05 X 10^ 


100 


2.80 X 10^ 


1.04 X 10^ 



over a Fermi-Dirac spectrum for given temperature. For 
our purpose, it is advantageous to have the total cross 
sections as a function of neutrino energy. In addition, 
charged-current cross sections for antineutrinos on tri- 
tons are needed. Therefore, we have computed the rele- 
vant cross sections using a model based on the random 
phase approximation (RPA) which has been successfully 
applied to the study of many neutrino-induced reactions 
(for example, see Ref. [33l.j34l|). Our approach follows 
the one described in Ref. f3S|, where we adopt an RPA 
that distinguishes between proton and neutron degrees 
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of freedom for the particle and hole states. The parent 
ground state is approximated by the lowest independent 
particle model state with the single-particle energies de- 
rived from an appropriate Woods-Saxon potential that 
reproduces the particle separation energies in the parent 
nucleus. The partial occupancy formalism as described in 
Ref. [35I is applied to the proton and neutron states and 
holes for ■^H and '^He, respectively. As residual interac- 
tion, we have used the Landau-Migdal force of Ref. [35[ . 
The calculation includes all multipole transitions with 
A ^ 4 and both parities, properly accounting for the 
dependence of the multipole operators on the momen- 
tum transfer [H, [s^l- The Gamow- Teller strength has 
been quenched by a factor 0.74 as in shell-model calcu- 
lations [3^. Our calculated cross sections are listed in 
Table |I] as a function of incident antineutrino energy E. 
To compare our results with the ab-initio calculations of 
Ref. 8], we have folded the cross sections with a Fermi- 
Dirac neutrino distribution and find agreement to better 
than a few percent in all cases. 

From the effective opacity of Eq. we determine the 
optical depth as: 



r(r) 



(3) 



The neutrinospheric radius Ri, is then defined as the po- 
sition where the optical depth reaches t{R^) = 2/3 |3g |. 
We assume that the neutrino distribution function at the 
neutrinosphere is represented by a Fermi-Dirac spectrum 
of temperature T and effective degeneracy parameter r/^: 



dE^ ' 2n^ (hc)^ exp[E/{kT^) 



1 



(4) 



The neutrino number luminosity [40] at the neutri- 
nosphere is then given by 



2c 



iknfRlfiR,)F2M 



(5) 



and a similar equation for antincutrinos. Here / is the 
flux factor needed to convert number density to number 
flux, with f{Rv) ~ 0.25 at the location of the neutri- 
nosphere [401 . The Fermi integral for relativistic particles 
is defined as Fn{vi) — dxx'^/[l + exp(a; — 77)], where 
77 should be regarded as a spectral parameter that is not 
necessarily related to the neutrino degeneracy. 

The above procedure allows for a completely indepen- 
dent determination of the electron neutrino and antineu- 
trino spheres, and consequently of their number luminosi- 
ties. However, the number luminosities are constrained 
by the requirement that the net flux of electron neutri- 
nos minus electron antincutrinos is positive or zero to 
warrant that the lepton content of the star does not in- 
crease. Simulations [4l[ show that the nascent neutron 
star deleptonizes only very gradually and therefore the 
number luminosities of neutrinos and antincutrinos are 
rather similar, allowing us to make the approximation 
that the net flux is equal to zero, Ln^u = Ln^p. This 



choice is consistent with our assumption of beta equilib- 
rium for fixing the composition in the outer layers of the 
neutron star (see Sect. Ill A)) , in which case the same num- 
bers of electron neutrinos and antincutrinos are created 
in this region. Since our calculations are mainly sensitive 
to the difference in the spectral properties of neutrinos 
and antincutrinos, we set ryp = for simplicity and obtain 
rji, from the condition of zero net flux. 

We use the following iterative procedure to deter- 
mine the Ve and Ve neutrinospheres. First, we assume 
some initial values for the temperatures of the neutrino 
and antineutrino distributions, and determine r]i, from 
the condition of zero net flux. Using these values we 
compute the radii, at which the optical depths become 
2/3. Then, we set the new and De temperatures to 
the local temperature at the corresponding radii and up- 
date 771/ from the zero net flux condition. We iterate this 
procedure until convergence. 

Since mainly antincutrinos are affected by the changes 
in composition, wc focus in Fig. [3] on the different con- 
tributions to the antineutrino absorption opacity (Kabs) 
and the effective opacity (Kcff), at four different times 
post bounce. The dashed lines represent our reference 
calculation, which considers only antineutrino absorption 
on protons, nucleon-nucleon bremsstrahlung and neglects 
the presence of light nuclei, corresponding to the treat- 
ment in Ref. [3] . For comparison, the solid lines are based 
on the virial EOS composition and the antineutrino inter- 
actions as discussed above (similar results are obtained 
using the NSE EOS). We find in Fig. [3] that the an- 
tineutrino (Kabs) in the region around the neutrinosphere 
is dominated by the contribution from deuterons and 
tritons, for which charged-current and neutral-current 
reactions contribute approximately in equal amounts. 
Although tritons become more abundant (see Fig. 
deuterons become more important for the opacities. 



III. NEUTRINO LUMINOSITIES, SPECTRA, 
AND NEUTRON EXCESS IN THE EJECTA 

In the following, we explore the influence of light nu- 
clei on the properties of the emitted neutrinos. To this 
end, we consider four different cases denoted A to D in 
Table |TT1 which use the different EOS discussed above. 
Our reference calculation is case A, which is performed 
with the same EOS (consisting of neutrons, protons, and 
a-particles in NSE) and the same neutrino reactions as in 
Ref. y\. When we change to the improved EOS including 
light nuclei for given density, temperature, and Y^, the 
different baryonic composition leads to a modified opti- 
cal depth and therefore to a shift of mainly the position 
of the neutrinosphere of electron antincutrinos. On the 
other hand, for given density and temperature in the neu- 
tron star atmosphere, the improved EOS imply different 
nucleon chemical potentials and therefore yield a differ- 
ent value of Ye when we impose the constraint of beta 
equilibrium for the stellar matter (see Sect. Ill A)) . This 
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FIG. 3: (Color online) Different contributions to the absorp- 
tion opacity (fCabss) and effective opacity (Keff) of electron an- 
tineutrinos as a function of radius for t = 2,5,7 and 10 s post 
bounce (from top to bottom) and for two different EOS: the 
one used in Ref. [7| (dashed lines) and the virial EOS (solid 
lines), based on the target abundances displayed in Fig. [5] 
In addition, the total absorption ("abs") and total effective 
opacity ("effective") are shown. The vertical lines mark the 
positions of the neutrinospheres of electron antineutrinos for 
the two different EOS. In the second panel from the top the 
vertical dashed and solid lines coincide. 



TABLE II: Different cases explored in Sect. |llll 



Case 


Ye determined from 


EOS and composition 


A 


beta equilibrium 


NSE (n, p, *He) 


B 


case A 


NSE (nucleons and nuclei) 


C 


beta equilibrium 


NSE (nucleons and nuclei) 


D 


beta equilibrium 


virial (n, p, A ^ 4 nuclei) 



0.25 r 



0.20 r 
0.15 : 

>- 

0.10 - 




0.05 r - 

n nn F - 
9.0 9.5 10.0 10.5 11.0 

rodlus [km] 

FIG. 4: Profile of the electron fraction Ye in the region around 
the neutrinosphere at f = 2s after core bounce. The dashed 
line corresponds to the EOS of Ref. (case A), the dot- 
ted line to the NSE EOS (case C), and the solid line to the 
virial EOS (case D). All profiles are obtained assuming beta 
equilibrium for the corresponding EOS. 

again influences the composition, optical depth, and neu- 
trinospheric positions. To quantify separately the impact 
due to the direct compostion change and the one associ- 
ated with an adjustment to a new beta equilibrium, we 
define an intermediate case B that uses the same den- 
sity, temperature, and profiles as case A but obtains 
the baryonic composition from the improved NSE EOS. 
Case C then uses the same EOS as case B but takes 
into account the adjustment to a new beta equilibrium. 
In case D, the baryonic composition as well as the beta 
equilibrium state are based on the virial EOS. 

Table IIIII shows the neutrinosphere radii and corre- 
sponding properties of electron neutrinos and antineu- 
trinos for the four different cases. The average energies, 
defined as 



and luminosities computed for our reference case A are 
in good agreement with the results of Ref. f7]. In or- 
der to compare the average energies and luminosities of 
Table IIIII with those commonly used in nucleosynthesis 
studies (which correspond to values measured at infin- 
ity), the gravitational redshift must be included. For our 
neutron star model, the redshift correction corresponds 
to a reduction by a factor of about 0.8. 

The differences in the various observables of Table IIIII 
caused by the presence of light nuclei can be understood 
from a detailed comparison of the four cases. While 
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TABLE III: Neutrinosphere radii Rv^,v^, neutrino spectral temperatures T;^^,^^, and average energies {^v^,v^), as well as number 
luminosities L„, spectral parameter 77^^, and wind electron fractions at four different times post bounce. 
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= 2s 










A 


10.01 


8.14 


25.64 


6.05 




0.72 


10.55 


6.34 


20.71 


0.514 


B 


9.977 


8.30 


26.16 


6.38 




0.79 


10.55 


6.34 


20.80 


0.507 


C 


10.00 


8.17 


25.73 


6.10 




0.73 


10.55 


6.35 


20.75 


0.513 


D 


9.979 


8.29 


26.12 


6.36 




0.77 


10.53 


6.37 


20.87 


0.509 












t 


= 5s 










A 


9.272 


7.17 


22.60 


3.55 




1.01 


9.821 


5.14 


17.10 


0.478 


B 


9.260 


7.24 


22.83 


3.65 




1.04 


9.819 


5.15 


17.16 


0.475 


C 


9.295 


7.04 


22.17 


3.37 




0.94 


9.814 


5.16 


17.07 


0.487 


D 


9.272 


7.17 


22.60 


3.55 




1.00 


9.813 


5.16 


17.15 


0.480 












t 


= 7s 










A 


9.107 


6.88 


21.69 


3.03 




1.15 


9.683 


4.73 


15.90 


0.462 


B 


9.095 


6.97 


21.95 


3.13 




1.19 


9.681 


4.74 


15.96 


0.458 


C 


9.139 


6.68 


21.04 


2.78 




1.04 


9.676 


4.75 


15.82 


0.475 


D 


9.134 


6.71 


21.14 


2.82 




1.05 


9.675 


4.75 


15.85 


0.473 












t ■- 


= 10s 










A 


9.041 


6.94 


21.86 


3.06 




1.49 


9.592 


4.37 


15.05 


0.431 


B 


9.039 


7.02 


22.12 


3.17 




1.53 


9.590 


4.37 


15.12 


0.427 


C 


9.063 


6.49 


20.44 


2.51 




1.23 


9.582 


4.39 


14.82 


0.456 


D 


9.065 


6.45 


20.32 


2.47 




1.20 


9.581 


4.39 


14.80 


0.458 



case A considers the neutron star atmosphere to be in 
beta equihbrium with only neutrons, protons, and alpha 
particles, case B uses the same Yf. profile as case A, but 
accounts for the presence of light nuclei. Their appear- 
ance happens at the expense of the number of free pro- 
tons, whose mass fraction is drastically reduced. This 
leads to a lower antineutrino opacity than in case A, 
even when the additional antineutrino interactions with 
light nuclei are fully included. Consequently, the neu- 
trinosphere of antineutrinos moves to a smaller radius 
where the temperature is larger so that antineutrinos are 
expected to be radiated with a higher mean energy. 

However, as discussed above, the EOS with light nu- 
clei also lead to a shift of the beta equilibrium conditions 
in the neutron star atmosphere. The balance of electron 
neutrino and antineutrino absorption and production re- 
actions is locally established for a higher value of Ye, thus 
compensating for the reduced abundance of free protons 
in the presence of light nuclei. This effect is shown in 
Fig.dlfor the Yg profile at t = 2s from the different EOS. 
The shift of the beta equilibrium is taken into account in 
cases C and D. As a consequence of the higher Fgi the 
abundances of protons, deuterons, and tritons are larger 
than in case B (although free protons are clearly reduced 
relative to their abundance in case A, see Fig. [5]). Com- 
pared to case B, the antineutrino opacity is therefore in- 
creased and the corresponding neutrinosphere is located 
at a larger radius, leading to a lower mean energy of the 
escaping antineutrinos. 

In contrast to the properties of electron antineutrinos, 
the spectra of electron neutrinos are only slightly affected 
by the improved EOS and by the adjustment to a new 



beta equilibrium. There are two reasons for this: The 
mass fraction of free neutrons dominates the composition 
and only slightly differs for all cases, and second, weak 
reactions with neutrons are responsible for most of the 
opacity of electron neutrinos. 

From this general discussion we conclude that the 
change of the baryonic composition due to light nuclei on 
the one hand and due to the neutrino-driven adjustment 
of matter to a new beta equilibrium on the other have 
effects in opposite directions for the position of the an- 
tineutrinosphere (and rather unimportant effects on the 
neutrinosphere). This makes it difficult to reliably pre- 
dict the change of this position in cases C and D com- 
pared to the reference case A. In fact, as Table Hill shows . 
the outcome of these competing effects can go either way. 
At early times {t = 2s after bounce) the matter in the 
vicinity of the neutrinosphere has large temperatures and 
hence a rather large fraction of free protons is present (see 
Fig. [2|) . Under such conditions protons dominate the an- 
tineutrino opacity in the region around the antineutri- 
nosphere (see Fig.[3l). In cases C and D the proton mass 
fraction is noticeably reduced, but the additional opac- 
ity contributions due to light nuclei cannot compensate 
for the reduction of the proton opacity. As a result the 
antineutrinosphere moves to slightly smaller radii result- 
ing in larger average energies for the emitted antineutri- 
nos. As the protoneutron star and its atmosphere cool, 
the surface density profile steepens (see Fig. [T]), and the 
neutrinospheres move to smaller radii within the same 
model. The matter in the region of the neutrinospheres 
becomes more neutron-rich and the proton abundance is 
lower. However, in cases C and D substantial amounts 
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of deuterons and tritons are present in this region over- 
compensating the reduction in the proton mass fraction 
and making the total antineutrino absorption opacities 
higher in cases C and D compared to case A. Conse- 
quently, the antineutrinosphere moves to larger radii re- 
sulting in smaller antineutrino average energies. 

We emphasize that at late times the antineutri- 
nosphere is located at densities p > lO^^g/cm'^, where 
nuclear interactions and many-body contributions affect 
the composition and neutrino cross sections. For exam- 
ple, for t = 7 and 10 s, the contribution from densities 
p > lO^^g/cm^ to the optical depth is 52% and 65%, 
respectively. In the present work, as well as in state- 
of-the-art studies of protoneutron star winds 0, such 
potentially important effects have been neglected. They 
should be considered in future work. 

In order to quantify the effects of the changing an- 
tineutrino energies on the nucleosynthesis conditions in 
the baryonic wind driven by neutrino energy deposition 
off the neutron star surface, we estimate the wind elec- 
tron fraction from the expression [ist 



— 



A,. 



(7) 



Here A^^„ and 



Al/eP 



are the neutrino absorption rate 
on neutrons and the antineutrino absorption rate on 
protons, respectively, which depend on the neutrino 
number luminosities Ln^^, the neutrino spectra, and 
the radial distance r from the neutron star: A^, = 
-^n,i/ (o'}/[47rr^/(r)], where (a) is the relevant cross sec- 
tion suitably averaged over the neutrino spectrum [1, 
and including weak magnetism corrections [29| . In mak- 
ing use of Eq. d?]), we assume that the ejected matter is 
initially composed only of neutrons and protons, so we 
neglect the presence of alpha particles and the so-called 
alpha-effect 4J|. Furthermore, we suppose that the mat- 
ter is exposed long enough to neutrino and antineutrino 
captures to achieve an equilibrium between neutrino and 
antineutrino absorptions, and that this happens at such 
large distances (and low temperatures) that electron and 
positron captures can be ignored. 

Using Eq. ([7]), we have calculated the wind electron 
fraction Y^ for the different cases. The corresponding 
results are given in Table IIIII Because of the increase 
of the average antineutrino energies, we find that matter 
is ejected in cases C and D at early times with slightly 
lower Y^ values than in case A. At later times {t > 5s 
after bounce) the mean antineutrino energies are smaller 
in cases C and D than in case A, and therefore YJ' is 
slightly higher than in our reference case. 

Once the ejected matter has reached larger radii and 
thus low enough temperatures, nuclei can form. During 
the early phases of nucleosynthesis mainly alpha parti- 
cles but also light nuclei are present in the composition. 
Neutrino interactions with these light elements can con- 
stitute an additional source of energy deposition in the 
wind. As noted by Qian and Woosley y such an ad- 
ditional source of energy could increase the entropy and 



reduce the dynamical time scale of the wind, and conse- 
quently facilitate the production of heavy nuclei via the 
r-process [ll]. Using one of the wind trajectories result- 
ing from the hydrodynamical simulations of Ref. ff\ , we 
have computed the contribution to the energy deposition 
rate arising from light nuclei (deuterons, tritons, and al- 
pha particles). Among the light nuclei, we found that 
the dominating contribution comes from neutrino inter- 
actions with alpha particles (where we have taken the 
cross sections from Gazit and Barnea jl^). However, our 
calculation showed that the additional energy deposition 
provided by light elements is too low, by more than an 
order of magnitude, to have an impact on the conditions 
for r-process nucleosynthesis. 



IV. CONCLUSIONS 

We have shown that the thermodynamical conditions 
in the outer layers of a protoneutron star favor the pres- 
ence of light nuclei, mainly deuterons and tritons, which 
are not accounted for by EOS currently used in core- 
collapse supernova simulations [isj and in studies of 
neutrino-driven supernova outflows Q. Using the pro- 
files of a hydrodynamical model for neutrino-driven su- 
pernova ejecta [7], we have estimated the effects of light 
nuclei on the emission of electron neutrinos and antineu- 
trinos. For this purpose we have compared the virial and 
NSE EOS, which include light nuclei, to a reference case 
composed of neutrons, protons, and alpha particles. 

The abundance of light nuclei can be studied in lab- 
oratory heavy-ion collisions. These experiments can re- 
produce the densities and temperatures near the neutri- 
nosphere. In the future, because the neutrinosphere is 
neutron rich, the abundances of light nuclei should be 
measured for more neutron-rich systems. This can be 
done with radioactive beams. 

The appearance of light nuclei has only a minor im- 
pact on the position of the electron neutrinosphere and 
consequently on the average energy of the radiated elec- 
tron neutrinos. However, the situation is different for 
electron antineutrinos. At early times when the pro- 
toneutron star is relatively hot and protons have mass 
fractions around 0.1 in the neutrinospheric region, the 
appearance of light nuclei reduces the antineutrino opac- 
ity. Therefore, antineutrinos escape from hotter layers 
in the protoneutron star with slightly larger average en- 
ergy. At later times the mass fraction of protons in the 
protoneutron star is greatly reduced, but light elements 
(in particular tritons) can have mass fractions that reach 
values even around 0.1. This makes light nuclei the major 
source of opacity for antineutrinos. Comparing with the 
EOS used in Ref. Q, we see that antineutrinos arc kept 
in thermal equilibrium with matter until larger radii in 
the protoneutron star, reducing the average energy of the 
emitted antineutrinos. For the latest time considered in 
the present study {t = 10s after bounce), the reduction 
could be as large as 1.5 MeV. 
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The changes in the antineutrino average energies can 
have consequences for the nucleosynthesis occurring in 
neutrino-driven winds. Such winds are a very interesting 
nucleosynthesis site. They allow for proton-rich ejecta 
during the first couple of second s l46l . |47| , where the re- 
cently suggested i/p-process [H, l49l| may take place, as 
well as neutron-rich ejecta at later times, which might 
provide the conditions for r-process nucleosynthesis [50j. 
During the early proton-rich phase, the changes in the 
antineutrino energies are minor and consequently the 
electron-to-baryon ratio in the wind remains the same 
or is slightly reduced. For the neutron-rich phase and in 
particular for the latest times considered in the present 
study, we find that can increase by as much as 0.025. 
This is a substantial change, and if everything else re- 
mains the same, such a change makes the occurrence of 
strong r-processing less likely. 

For more detailed studies, reliable estimates of the neu- 
trino cross sections with light nuclei are desirable. The 
neutrino-deuteron cross sections ^ used here are prob- 
ably sufficiently accurate at the neutrino energies con- 
sidered in this work. For charged-current and neutral- 
current cross sections of antineutrino reactions with tri- 
tons, we have presented results based on a relatively sim- 
ple RPA approach. The neutral-current cross sections 
were found to agree very well with the recent ab-initio re- 
sults of Ref. Q . Similar calculations for charged-current 
cross sections would be useful. 

Our estimates show that future simulations of 
neutrino-driven supernova outflows should take into ac- 
count light elements in the baryonic composition of 
the stellar medium and the corresponding cross sec- 



tions of neutrino interactions, especially those of elec- 
tron antineutrinos. This is particularly important to 
fully quantify the consequences for the properties of the 
neutrino emission from forming neutron stars and for 
the neutrino-generated nucleosynthesis conditions in the 
baryonic mass that is lost from such stars. In addition 
to including light nuclei in the EOS, our results show 
that, for late times, an important future problem is un- 
derstanding the properties of and neutrino interactions 
with nucleonic matter at densities above 10"'^^ g/cm'^. 
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